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ABSTRACT 

HD 133880 is a rapidly rotating chemically peculiar B-type (Bp) star (u sin i ~ 103 km s -1 ) 
and is host to one of the strongest magnetic fields of any Ap/Bp star. A member of the 
Upper Centaurus Lupus association, it is a star with a well-determined age of 16Myr. 12 new 
spectra, four of which are polarimetric, obtained from the FEROS, ESPaDOnS and HARPS 
instruments, provide sufficient material from which to re-evaluate the magnetic field and obtain 
a first approximation to the atmospheric abundance distributions of He, O, Mg, Si, Ti, Cr, Fe, 
Ni, Pr and Nd. 

An abundance analysis was carried out using zeeman, a program which synthesizes spectral 
line profiles for stars with permeating magnetic fields. The magnetic field structure was char- 
acterized by a colinear multipole expansion from the observed variations of the longitudinal 
and surface fields with rotational phase. Both magnetic hemispheres are clearly visible during 
the stellar rotation, and thus a three-ring abundance distribution model encompassing both 
magnetic poles and magnetic equator with equal spans in colatitude was adopted. 

Using the new magnetic field measurements and optical photometry together with previously 
published data, we refine the period of HD 1 33880 to P = 0.877 476 ± 0.000 009 d. Our simple 
axisymmetric magnetic field model is based on a predominantly quadrupolar component that 
roughly describes the field variations. Using spectrum synthesis, we derived mean abundances 
for O, Mg, Si, Ti, Cr, Fe and Pr. All elements, except Mg, are overabundant compared to 
the Sun. Mg appears to be approximately uniform over the stellar surface, while all other 
elements are more abundant in the negative magnetic hemisphere than in the positive magnetic 
hemisphere. In contrast to most Ap/Bp stars which show an underabundance in O, in HD 
133880 this element is clearly overabundant compared to the solar abundance ratio. 

In studying the Ha and Paschen lines in the optical spectra, we could not unambiguously 
detect information about the magnetosphere of HD 133880. However, radio emission data 
at both 3 and 6 cm suggest that the magnetospheric plasma is held in rigid rotation with the 
star by the magnetic field and further supported against collapse by the rapid rotation. Subtle 
differences in the shapes of the optically thick radio light curves at 3 and 6 cm suggest that the 
large-scale magnetic field is not fully axisymmetric at large distances from the star. 
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♦Based in part on observations made with the European Southern Observatory (ESO) telescopes under ESO programmes 082.D-0061(A), 083.D-0034(A), 
085.D-0296(A) and 086.D-0449(A), obtained from the ESO/ST-ECF Science Archive Facility. It is also based in part on observations carried out at the 
Canada-Francc-Hawaii Telescope (CFHT) which is operated by the National Research Council of Canada, the Institut National des Science de TUnivers of 
the Centre National de la Recherche Scientifique of France and the University of Hawaii. 
fE-mail: jbaile33<?uwo.ca 


© 2012 The Authors 

Monthly Notices of the Royal Astronomical Society © 2012 RAS 


HD 133880 329 


1 INTRODUCTION 

HD 133880 (HR 5624) is a rapidly rotating late B-type chemi- 
cally peculiar (Bp) star with the Si A4200 peculiarity. It exhibits 
an unusual magnetic field. The mean line-of-sight magnetic field 
(B z ) is very' strong and observed to vary from about -4 to +2kG 
(Landstreet 1990). Unlike most chemically peculiar stars, HD 
133880 has a field that is predominantly quadrupolar as opposed 
to dipolar. HD 133880 is a photometric variable with variations of 
the order of 0. 15 mag in the U band, which may be the result of the 
large magnetic field (Waelkens 1985). The magnetic field serves to 
create a ‘patchy’ distribution of elements on the surface of the star. 
This non-uniform abundance distribution, and its associated line 
blocking and back warming, may explain the photometric variations 
(Ryabchikova 1991). Specifically, several elements are distributed 
non-uniformly over the stellar surface in a non-axisymmetric pattern 
about the rotation axis. As the star rotates, the observed magnetic 
field strength (fi z ) varies and over/underabundances of elements are 
detected. The most striking anomalous abundances in Bp stars are 
often found for Cr and rare earth elements which can, respectively, 
be as much as 10 2 and 10 4 " 5 overabundant compared to the Sun 
(Ryabchikova 1991). The resultant variability is explained by the 
oblique rotator model: the rotation and magnetic field axes are at 
angles i and fi to the line of sight and rotation axis, respectively. 

Landstreet ( 1 990) used longitudinal magnetic field measurements 
obtained from H p in conjunction with photometric data reported by 
Waelkens (1985) to deduce a period of P — 0.877 485 ± 0.00002 d 
for HD 133880. Lim, Drake & Linsky (1996) did a preliminary 
analysis of the rotational modulation of 3- and 6-cm radio emission 
from HD 133880. The radio emission variation had both broad and 
narrow peaks that correlated with the maximum extrema, respec- 
tively, of the dipolar, £ d , and quadrupolar, B qy contributions to the 
magnetic field reported by Landstreet (1990). 

Landstreet et al. (2007) confirmed HD 1 33880 as a member of the 
Upper Centaurus Lupus association and thus it is a star with a well- 
known age of log t = 7.20 ± O.lOyrs. HD 133880 is a very young 
star, having completed only about 5 per cent of its main-sequence 
lifetime (Landstreet et al. 2007, 2008). Their preliminary analysis of 
many physical characteristics concluded that = 1 2 000 ± 500 K, 

logZTLQ = 2.10 ± 0.1 and M/Mq — 3.20 ±0.15. [Our analysis 
(see below) has shown that the T& is likely closer to 13 000 K.) 
Table 1 summarizes the physical properties of HD 133880. 

There are few hot chemically peculiar stars with broad spectral 
features for which abundance analyses have been carried out. Be- 
cause of the very complex, strong and unusual magnetic field, HD 
133880 is a worthwhile candidate to study in detail. It is also in- 
teresting from the fact that it is a star with a well-known age. HD 
133880 is one in a large sample of cluster Ap/Bp stars first compiled 
by Bagnulo et al. (2006) and Landstreet et al. (2007, 2008) to study 
the time evolution of magnetic fields in Ap stars. The current study 
is part of an effort to understand empirically how chemical abun- 
dances evolve with time and stellar age in the magnetic chemically 
peculiar stars. Recently, Bailey et al. (201 1) characterized one star 
in this sample, the Ap star HD 318107 (NGC 6405 77). HD 133880 
is a second star of great interest from that sample. 

This paper will discuss the efforts to characterize the mag- 
netic field and abundance distributions of several elements for HD 
133880. The following section discusses the polarimetric and spec- 
troscopic observations used as well as the deduced longitudinal 
fields; Section 3 discusses improvement of the rotational period; 
Section 4 describes the determination of physical parameters; Sec- 
tion 5 addresses the observed line profile variations; Sections 6 and 


Ifcble 1. Summary of stellar, wind, magnetic and magnetosphcric properties 
of HD 133880. 


Spectral type 

B8IVp Si X4200 

SIMBAD and BSC 

TeSf (K) 

13000 ± 600 

This paper 

log 8 (cgs) 

4.34 ±0.16 

This paper 

R* (Ro) 

2.01 ± 0.32 

This paper 

vsin i (kms -1 ) 

103 ± 10 

This paper 

P(d) 

0.877 476 ±0.000 009 

This paper 

log (£*/Lq) 

2.02 ±0.10 

This paper 

(Mo) 

3.20 ± 0.15 

Landstreet et al. (2007) 

log t (Myr) 

7.20 ±0.10 

Landstreet et al. (2007) 

Rd(G) 

-9600 ± 1000 

This paper 

R q (G) 

-23000 ± 1000 

This paper 

Boct (G) 

1900 ± 1000 

This paper 

to 

55 ± 10 

This paper 

pn 

78 ± 10 

This paper 

v* 

~10 7 

This paper 

w 

0.2 

This paper 

Tfpin 

11 Myr 

This paper 


7 discuss the magnetic field and abundance models; Section 8 re- 
ports the abundances obtained; Section 9 discusses the Hof and radio 
magnetosphere and Section 10 summarizes the presented work. 

2 OBSERVATIONS 
2.1 Polarized spectra 

High-resolution spectropolari metric (Stokes / and V) observations 
of HD 1 33880 were collected with ESPaDOnS at the 3.6-m Canada- 
France-Hawaii Telescope (CFHT) and HARPSpol on the Euro- 
pean Southern Observatory (ESO) 3.6-m telescope. ESPaDOnS and 
HARPSpol are both high-resolution fibre-fed spectropolarimeters. 
Three ESPaDOnS observations were obtained on 201 0 August 3 and 
on 201 1 July 14 and 15 [the latter in the context of the Magnetism 
in Massive Stars Large Programme (MiMeS)], and one HARPSpol 
spectrum was obtained on 20 10 May 3. Each ESPaDOnS spectropo- 
larimetric sequence consisted of four individual subexposures taken 
in different retarder configurations, whereas the HARPSpol spec- 
trum results from the combination of four independent spectra, each 
comprised of four subexposures, obtained over the course of about 
9 min. The ESPaDOnS measurements cover a spectral range from 
370 to 1050 nm at a spectral resolution of R ~ 65000, while the 
HARPSpol observation covers 370-690 nm at a spectral resolution 
of/? - 100000. 

From each set of four subexposures we derived Stokes / and 
Stokes V spectra following the double-ratio procedure described by 
Donati et al. (1997), ensuring in particular that all spurious signa- 
tures are removed at first order. Null polarization spectra (labelled 
N) were calculated by combining the four subexposures in such 
a way that polarization cancels out, allowing us to verify that no 
spurious signals are present in the data (see Donati et al. 1997, for 
more details on the definition of AO- All ESPaDOnS frames were 
processed using the CFHT’s Upena pipeline, feeding the automated 
reduction package Libre-Esprit (Donati et al. 1997). The HARP- 
Spol spectrum was reduced using the reduce code (Piskunov & 
Valenti 2002). The details of this reduction can be found in recent 
papers by Makaganiuk et al. (2011) and Kochukhov et al. (201 lb). 
The peak signal-to-noise ratios (S/N values) per 1.8 km s _1 velocity 
bin in the reduced ESPaDOnS and HARPSpol spectra range from 
449 to 697. 
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liable 2. Log of ESPaDOnS and HARPSpol observations of HD 133880. Listed are the instrument used, the 
Heliocentric Julian Date of the mid-point of the observation, total exposure time, the peak S/N per 1.8kms -1 
velocity bin, the phase of the observation (according to equation 2), the evaluation of the detection level of a Stokes 
V Zeeman signature (DD = definite detection, MD = marginal detection and ND = no detection) and the derived 
longitudinal field and longitudinal field detection significance z from both V and N. In no case is any marginal or 
definite detection obtained in the N profiles. 


Instrument 

HJD 

*exp 

(s) 

S/N 

pixel -1 

Phase 

Detect 

V 

Bt ± ob 
(G) 

z 

N 

Bt ± cr B 
(G) 

z 

ESPaDOnS 

2455411.732 

720 

629 

0.640 

DD 

2029 ± 38 

53.4 

-51 ±31 

1.6 

ESPaDOnS 

2455756.795 

720 

697 

0.878 

DD 

-2103 ± 38 

55.3 

43 ±23 

1.9 

ESPaDOnS 

2455757.877 

720 

573 

0.111 

DD 

-3671 ± 43 

85.4 

-23 ± 26 

0.88 

HARPSpol 

2455319.726 

520 

449 

0.787 

DD 

650 ± 62 

10.5 

-15 ±58 

0.26 


Table 3. Log of unpolarized spectra. Listed arc the in- 
strument used, the Heliocentric Julian Date, the peak S/N 
per 1.8kms -1 velocity bin, the exposure time and the 
phase computed using the ephemeris of equation (2). 


Instrument 

HJD 

S/N 

fexp (s) 

Phase 

ESPaDOnS 

2455412.727 

180 

386 

0.774 

ESPaDOnS 

2455414.728 

180 

394 

0.054 

FEROS 

2454867.846 

130 

143 

0.809 

FEROS 

2454868.818 

130 

103 

0.917 

FEROS 

2455031.608 

250 

122 

0.438 

FEROS 

245 5599.820 

240 

141 

0.990 

FEROS 

245 5600.828 

240 

151 

0.139 

FEROS 

245 5601.868 

300 

121 

0.325 


The leg of spectropolarimetric observations is presented in 
Table 2. 


2.2 Unpolarized spectra 

High-resolution spectroscopic observations (Stokes 7) were ob- 
tained from ESPaDOnS at CFHT and ESO’s FEROS instrument 
at the 2.2- m telescope located at La Silla. FEROS is an echelle 
spectrograph that has a wavelength coverage from 350 to 920 nm 
with a spectral resolution of R ~~ 48 000. FEROS is fed by two 
fibres (object and sky fibres), which in principle can be used to 
subtract the contribution of the sky from the object spectrum. A 
more detailed description is given by Kaufer et al. (1999). The two 
ESPaDOnS observations were taken on 2010 August 4 and 6, and 
the six FEROS spectra were collected on 2009 February 5, 6 and 
July 19, as well as on 201 1 February 7, 8 and 9. 

The ESPaDOnS data were reduced using Libre-Esprit in a man- 
ner similar to that described above. The FEROS data were reduced 
using the data-reduction software implemented under mjdas (Kaufer 
et al. 1999). The log of unpolarized spectra is presented in Table 3. 


This process allows the computation of single averaged Stokes I 
and V profiles with much higher S/N values than those of the in- 
dividual lines, dramatically improving the detectability of Zeeman 
signatures due to stellar magnetic fields. 

The atomic data were taken from the Vienna Atomic Line 
Database (VALD; Piskunov et al. 1995; Ryabchikova et al. 1997; 
Kupka et al. 1999, 2000). The VALD line list used in this analysis 
is characterized by the stellar effective temperature T c fr and surface 
gravity logg reported in Table 1, abundances as derived from the 
detailed abundance analysis (see Section 8) and an unbroadened line 
depth threshold equal to 1 per cent of the continuum intensity 7 C , 
yielding approximately 10000 lines in the ESPaDOnS spectral win- 
dow and 8000 lines in the HARPSpol window. The line depths were 
then interactively adjusted so that the line depths predicted by the 
LSD convolution model matched the observed depth of the spectral 
absorption lines as well as possible, thus improving the fit between 
the spectrum and the LSD model. Because Stokes V profiles are de- 
tected in many individual spectral lines of HD 1 33880, we verified 
(by examining simultaneously the Stokes V and 7 profiles) that this 
adjustment procedure also resulted in a better agreement between 
the observed and computed Stokes V spectrum. 

The mean longitudinal magnetic field was measured by comput- 
ing the first-order moment of the Stokes V LSD profile within the 
line according to the expression 


<S Z ) = -2.14 x 10“ 


/i>V(u)dv 

*Scffc/[/ c - /(u)]du' 


( 1 ) 


where {B z ) is in G, g eff is the effective Landg factor, A. is the mean 
wavelength of the LSD line in A and v is the velocity within the 
profile measured relative to the centre of gravity (Donati et al. 1997; 
Wade et al. 2000). The uncertainties associated with (B r ) were deter- 
mined by propagating the formal (photon statistical) uncertainties 
of each pixel through equation (1). The integration ranges employed 
in the evaluation of equation (1) associated with each Stokes profile 
were selected individually through visual inspection so as to include 
the entire span of the (highly variable) Stokes I profile. 


23 Longitudinal magnetic field measurements 

The circular polarization observations were analysed using the mul- 
tiline analysis technique least squares deconvolution (LSD; Donati 
et al. 1997). LSD combines the information from essentially all 
metallic and He lines in the spectrum by means of the assumption 
that the spectrum can be reproduced by the convolution of a single 
‘mean’ line profile (the LSD profile) with an underlying spectrum 
of unbroaiened atomic lines of specified line depth, Landd factor 
and wavelength (the line mask, as described by Wade et al. 2000). 


2.4 ATCA radio continuum measurements 

HD 133880 was observed with the Australia Telescope Compact 
Array (ATCA) at 6- and 3.5-cm wavelength simultaneously on 
1995 February 12, 14 and 16. On each day, the observation spanned 
~10h, during which time we observed HD 133880 continuously 
apart from short (~~3 min) scans of a secondary calibrator at regular 
(^20 min) intervals. In this way, we were able to attain full cover- 
age in rotation phase for the star, with repeated coverage (on two 
separate days) for just over half a rotation phase to check where 
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any observed variations were repeatable (as indeed they were, as 
shown in fig. 2 of Lim et al. 1996). The primary (absolute flux 
density) calibrator used was PKS 1934—638, and the secondary 
(complex gain) calibrator PKS 1458—391. We reduced the data in 
the standard fashion using miriad, being careful to omit data with 
poor phase coherence (typically at the beginning and end of each 
day when the source was low in the sky). 

3 ROTATION PERIOD 

A period of 0.877 46 ± 0.000 0 1 d was reported by Waelkens (1985) 
based on Geneva photometric observations. Landstreet (1990) rec- 
ommenced an adjustment of the period to 0.877 485 ± 0.000 020 d 
to bring the magnetic measurements reported by Borra & Landstreet 
(1975) into agreement with those that he had acquired in 1987 and 
1988. 

When we phase the accumulated magnetic measurements of 
HD 133880 [i.e. the three measurements of Borra & Landstreet 
(1975), the 12 measurements of Landstreet (1990) and the four new 
measurements presented here] according to the Landstreet (1990) 
ephemeris, we observe a systematic offset of ~0. 1 cycles between 
the new measurements and those published previously. Such an off- 
set could be interpreted as the consequence of a small error in the 
assumed period. However, given the presence of important abun- 
dance nen-uniformities on the surface of HD 1 33880, it is likely that 
the shape and possibly the amplitude of the longitudinal field vari- 
ation derived from metallic lines (i.e. our new LSD measurements) 
and those derived from hydrogen lines (i.e. the previously published 
measurements) may differ. Given the sparse phase coverage of our 
new measurements, it is therefore ambiguous whether this apparent 
offset is truly a consequence of a period error. 

In an attempt to resolve this uncertainty, we examined the pho- 
tometric measurements obtained by the Hipparcos mission (ESA, 
1997). 130 7/p-band measurements of HD 133880 are reported 
in the 1997 edition of the catalogue. 1 The periodogram of those 
measurements in the range 0.5-1. 5 d shows a clear, unique peak at 
0.877 479 ± 0.000 030 d. [The uncertainty on the period, which cor- 
responds to 1<t confidence, is computed assuming Gaussian statis- 
tics by calculating the variation of the reduced x 2 of a sinusoidal 
fit to the phased measurements in the vicinity of the best-fitting 
period (e.g. Press, Flannery & Teukolsky 1986)1. A similar analy- 
sis applied to 37 V-band measurements of Waelkens (1985) yields 
0.877 459^ ;^o 6 d. Both of these periods are consistent within the 
(lcr) error bars. To improve the precision of the period, we com- 
bined the two data sets. First, we used the transformation reported 
by Harmanec (1998) to convert the H p measurements to Johnson 
V-band. We then performed the same period analysis on the com- 
bined data set (spanning nearly 15 years), obtaining a period of 
0.877 4731 db 0.000 0060 d (see Fig. 1). 

Therefore the photometric data are indicative of a period of 
0.8774731 d, intermediate between those of Waelkens (1985) and 
Landstreet (1990). 

We note that if we phase the magnetic data with the above photo- 
metric period, a slightly longer period of 0.877 476 d is necessary to 
bring the new magnetic measurements into better agreement with 
the older measurements. Finally, this period also removes the phase 
offset between the observed radio flux variation noted by Lim et al. 
(1996) and the magnetic extrema. 


1 vizier.u-strasbg.fr/viz-bin/VizieR-S?HIP%2074066) 
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Figure 1. Phased V-band (Waelkens 1985; filled squares) and // p - band 
(ESA, 1997; open diamonds) photometry. The Hipparcos photometry has 
been scaled to the Johnson V system using the transformation of Harmanec 
(1998). 


Based on these results, we adopt the following ephemeris for HD 
133880: 

JDhct"- = 244 5472.000( 10) + 0.877 476(9)£, (2) 

where we have adopted a zero-point corresponding to minimum 
photometric brightness (which occurs simultaneously with the min- 
imum of the longitudinal magnetic field). We determined this (3or) 
uncertainty by examining the x 2 of the fit to the photometric data 
about this best-fitting value. 

4 DETERMINATION OF PHYSICAL 
PARAMETERS 

4.1 Effective temperature 

The enhanced metal abundances and the effects of the magnetic 
field have to be taken into account when determining the effective 
temperatures of Ap stars (Stepien & Dominiczak 1989; Hauck & 
Kunzli 1996). As a result, effective temperatures determined us- 
ing Geneva and/or Strbmgren uvbyP photometry using normal star 
calibrations must be corrected to the Ap stars' temperature scale. 
However, Hauck & North (1993) suggest that stars that are either 
He weak or He strong do not require that their temperatures be 
corrected. 

HD 133880 has available both Geneva and Stromgren uvbyP 
photometry for the determination of effective temperature. We have 
used both in this analysis. For the Geneva photometry we have 
used the fortran program described by Kunzli et al. (1997), and 
discussed by Landstreet et al. (2007), to obtain an estimate of the ef- 
fective temperature. For the Stromgren uvbyp photometry we used 
the fortran code ‘uvbybbtanew-ap’ which is a modified version 
of ‘uvbybetanew’ (Moon & Dworetsky 1985; Napiwotzki, Schoen- 
berner & Wenske 1993) that corrects the effective temperature to 
the Ap temperature scale. However, in this analysis no correction 
was applied to the temperature because HD 133880 is classified as 
a He- weak star (Hauck & North 1993), which is confirmed in this 
paper (see Section 8). 

We found T eff = 13 300 and 1 2 700 K using uvbyp and Geneva 
photometries, respectively. The scatter of these measurements about 
their mean is of the order of 300 K. Given the intrinsic uncertainties 
in determining the effective temperature for a magnetic Ap star 
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in addition to the global uncertainties in the applied methods, we 
optimistically adopt an uncertainty of approximately ±600 K. Both 
the Geneva and uvbyfi photometries are rotationally averaged values 
and their variations are included within the quoted uncertainties. 
We use the average of the computed temperatures, adopting an 
effective temperature of T& = 13000 ± 600 K for HD 133880. If 
the Ap correction to the temperature is applied, a value of T cff = 
1 2 440 K is found, which is within the uncertainty of our adopted 
value. 

Netopil et al. (2008) performed an analysis on the temperature 
calibration of chemically peculiar stars. Included in their analysis 
was HD 133880. Using photometry, their calibrated temperature 
corrections suggest that the T cff is 11930 ± 210 K. This agrees, 
within uncertainties, to our derived T cff if we corrected our derived 
temperature to the Ap temperature scale. In our analysis, we favour 
the higher T & based on work by Hauck & North (1993) and the 
better quality fits to H/3 using the model fits of Kurucz (1979) with 
the higher of these two temperatures. We do, however, compare 
the derived abundances of all elements using the 7 cff reported by 
Netopil et al. (2008) of 12 000 K to our adopted value of 13 000 K 
(see Section 8). 

4.2 Luminosity 

In addition to T & and the V magnitude of HD 133880 (see 
Table 1), we also require a bolometric correction in order to cal- 
culate the stellar luminosity. Landstreet et al. (2007) derived a set 
of bolometric corrections as a function of T eff that can be applied 
to Ap stars. Applying this correction (see equation 1 in Landstreet 
et al. 20C7) results in log ULq = 2.02 ±0.10. Here, the uncer- 
tainty was estimated in the same manner as described in Section 1 
and by Landstreet et al. (2007). Note that we used E(B —V) = 0.00 
to compute the luminosity, which is confirmed by looking at the 
colour-colour diagram of (C7 — B) versus (B — V). 


4.3 Other parameters 

Based on values of T& and luminosity determined in Sections 4.1 
and 4.2, the stellar radius is found to be R = 2,01 ± 0.32 Rq. 
Landstreet et al. (2007) report a mass for HD 133880 of M = 
3.20 ± 0.15 Mo (consistent with the values of T t and luminosity 
in this paper), which can be used with the derived radius to estimate 
the surface gravity, log g . We find log g — 4.34 ± 0. 1 6, consistent 
with young age. 

From our adopted period of P = 0.877 476 ± 0.000 009 d (Sec- 
tion 3) and v sin i = 103 ± 10 km s -1 (see below), we can estimate 
the inclination angle of the rotation axis to the line of sight, /, using 
the equation 


where v sin / is in km s" 1 and R is in Rq, to obtain a value of i = 
63 ± 18°. All the physical parameters derived are summarized in 
Table 1 and are in agreement with those reported by Landstreet et al. 
(2007). 

5 LINE PROFILE VARIATIONS 

To characterize the line profile variability of HD 133880, we mea- 
sured the equivalent width variations for a number of spectral lines. 
The spectral lines were renormalized to the surrounding continuum 
regions and the equivalent widths were computed by numerically 
integrating over the line profile. A single uncertainty value was 
estimated for each pixel from the rms scatter in the continuum re- 
gions surrounding the line profile, and the resulting uncertainties in 
the equivalent width measurements were found by adding the pixel 
uncertainties in quadrature. 

In Fig. 2 we show equivalent width measurements phased to the 
ephemeris given in equation (2) for select lines of titanium, iron. 



0.40 - 


£ 

w 

S 0.35 

lO 

*0) 0.30 


0.25 

0.65 

E 

0.60 

m 

° 0.55 


CO 


0.50 


ih 

Ih : 



. 

I 

i \ 

f 

| • 

: { 

i 

f 

f 

"fit 

t 

: 1 

* 


-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 
Phase (P= 0.877476 d) 


-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 
Phase (P= 0.877476 d) 


Figure 2. Phased equivalent widths measurements for the indicated spectral lines measured from the ESPaDOnS (black squares), FEROS (red triangles) and 
HARPSpo'i (blue circles) data sets. 


© 2012 The Authors, MNRAS 423, 328-343 
Monthly Notices of the Royal Astronomical Society © 2012 RAS 


HD 133880 333 


1-20 

Till 4533 

1.20 

Fell 4603 

* 

1.00 

0.80 

§fgg||| 

1.00 

O.BO 


0.80 

■ : 

0.60 

'm^'iiiiI ' liiir 

| 0/40 

£ 

0.20 


| 040 
0.20 


0.00 


0.00 


-0,20 


-0.20 



4S3Z 4933 4934 4339 4338 433? 
wavelength (angstrom*) 

Crn 4580/4592 


4382 4393 4604 4936 4589 

wavelength (angstroms) 

Sill 5056 



Figure 3. Observed metallic line profiles for the indicated lines. Shown 
are the continuum-normalized spectra, displayed in such a way that the 
continuum of each spectrum is plotted at a position on the vertical axis that 
corresponds to the phase of the observation. 


chromium and silicon. The equivalent width curves are all simi- 
lar and show strong, approximately sinusoidal variations consistent 
with the rotation period. For each element, the equivalent width 
reaches a minimum at about phase 0.5. The largest variation is seen 
in the Cm 4588 A line, with the equivalent width measurements al- 
most doubling in value. Similar equivalent width changes are found 
in the Tin 4533 A and Fe n 4583 A lines, but a much smaller ampli- 
tude of the variation is measured in the Si ii 5056 A doublet. Clear 
variations are also seen in other spectral lines of these dements. 

The equivalent width variations for each of the elements plotted 
in Fig. 2 appear to be the result of absorption features travelling 
through the line profile. The net result is a variation in the total 
depth, as well as shape, of the line, as demonstrated by the profiles 
shown in Fig. 3. The Si line plotted in Fig. 3 has the simplest 
variations, with what appears to be a single feature travelling from 
negative to positive velocities. The variations in the other lines are 
phased similarly, but the line profiles are generally more complex. 


rotation axis, /J, with a value of 90' (i = p = 90°). The dipolar and 
quadrupolar components were fi d = -81 25 G and B q = — 10 900 G, 
respectively. (Note that we have corrected negative sign errors in 
both components for the derived field reported by Landstreet 1990.) 
From this model, it was established that the predominant magnetic 
field distribution over the stellar surface is quadrupolar with (B z ) 
variations that are clearly not sinusoidal. Further, since a change 
in the sign of the longitudinal field is observed, i + p must be 
greater than 90°. However, Landstreet’s adopted geometry is barely 
consistent with our calculated inclination of i = 63 = ± 18 ° (see 
Section 4). Also, Landstreet ( 1 990) had no information about the 
variations of the surface magnetic field modulus, (fl), making it 
difficult to constrain the field further. Landstreet (1990) pointed out 
that there is not a unique model for the observed magnetic curve of 
HD 133880 and suggested an alternate geometry of i = 80° and p = 
55 c that produces a similar curve (note that the values of i and P can 
be interchanged without affecting the magnetic curve). This second 
geometry is therefore in much better agreement with our derived 
value of i. Unfortunately, Landstreet (1990) did not provide the 
values of the dipole and quadrupole components for this alternate 
geometry. We therefore search for a new magnetic field model that 
agrees well with our derived value for i of 63° ± 18°. 


6.1 Surface magnetic field modulus measurements 

To further constrain the magnetic field, information about the 
surface magnetic field modulus {B) is required. To obtain this, 
we cannot rely on Zeeman splitting because rotational broaden- 
ing dominates the spectrum, smearing out any observable Zeeman 
components. Although dominated by rotation, the magnetic field 
contributes a measurable fraction to the width of a (Stokes I) spec- 
tral line. Preston (1971) outlines a method to obtain mean surface 
field measurements from rotationally broadened spectral lines. By 
comparing the widths of two spectral lines of a given element that 
have a large and small mean Lande factor (i.e. lines that are strongly 
effected by the local magnetic field to ones that are not), we can 
obtain a meaningful measurement of the magnetic field modulus. 
Preston (1971) defines a parameter K : 


K = 


/ wg-wj \ 1/2 

^(XV) L -(XV)J * 


(4) 


where u;l and ws are the averages of the measured line widths, 
in cm, of lines with a large and small Lande factor, respectively, 
X is the wavelength in cm and z is the mean Lande factor where 
the subscripts L and S refer to the lines with large and small Lande 
factors, respectively. Preston then finds the relationship between the 
surface field and K to be 


B s (kG) = 0.5 4- 7.9 K. (5) 

The challenge with HD 133880 was to find a pair of spectral lines 
that were not strongly blended with other lines from which to mea- 
sure line widths. Table 4 presents two sets of lines for Cr ii and Fe ii 


6 MAGNETIC FIELD MODEL 

Before any meaningful abundance analysis can be done, we need 
to establish an appropriate magnetic field model. Landstreet (1990) 
derived a best-fitting magnetic field model based on longitudinal 
field measurements of Hp for HD 133880. He adopted a magnetic 
field geometry with the inclination angle of the line of sight to the 
rotation axis, /, equal to the angle of the magnetic field axis to the 


Table 4. Lines used to measure {B) . The atomic data 
were gathered from the VALD data base. 



Small 

z 

Large z 


Element 

*• (A) 

z 

mA) 

Z 

Fen 

4508.288 

0.50 

4520.224 

1.34 

Crn 

4284.188 

0.52 

4261.913 

1.080 
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Table 5. ( B ) measurements from lines of Cm and 
Fen for spectra from ESPaDOnS and HARPSpol. 
Phases are computed from equation (2) and uncer- 
tainties estimated to be ±2.5 kG. 




<B> (kG) 

Instrument 

Phase 

Fe ii 

Crn 

ESPaDOnS 

0.054 

19.4 

17.2 

ESPaDOnS 

0.111 

21.1 

22.9 

ESPaDOnS 

0.640 

10.7 

11.9 

ESPaDOnS 

0.774 

10.9 

12.8 

HARPSpol 

0.787 

11.4 

12.3 

ESPaDOnS 

0.878 

16.9 

14.8 


from which line widths could be measured. Although good candi- 
dates, these spectral lines are still quite blended with adjacent lines 
in the spectrum. To combat this problem, we measured line widths 
at the continuum by repeatedly fitting Gaussians to the spectral line 
using iraf’s splot function. The average width for each feature is 
taken in our calculation. We were also limited by resolving power 
and S/N in ourFEROS spectra. Specifically, the combination of the 
lower resolution of the FEROS instrument and lower S/N values in 
the data meant that we were unable to unambiguously measure the 
spectral line widths. As a result, only the ESPaDOnS and HARP- 
Spol spectra could be utilized. Results are summarized in Table 5 
and illustrated in the top panel of Fig. 4. It is encouraging that the 
magnetic field values measured from both Fen and Cm agree well 
with one another, suggesting the variations are indeed real. Due to 
the high blending in the spectral lines and the ambiguity in measur- 
ing the line widths, we estimate the uncertainties to be ~-±2500G, 
though this value may be optimistic. Although this method is not as 
precise as measuring the splitting of spectral lines due to the Zee- 
man effect, it does provide meaningful constraints on the magnetic 
field. It is clear from these measurements that the surface field varies 
from about 20 kG at the negative magnetic pole to 10 kG near the 
positive magnetic pole (see below). The derived field is consistent 
with a non-sinusoidal curve (though no definitive conclusions 
can be made given the limited number of data points and size of 
the error bars), much in the same manner to the Iine-of-sight mag- 
netic field, but in the opposite sense: the maximum of the longitudi- 
nal field curve corresponds to the minimum of the surface integrated 
field curve. More measurements of {£) are required to constrain the 
field further, which requires more high S/N observations of HD 
133880. 

6.2 Magnetic field geometry 

The fldsrch.f program was used to obtain the magnetic field geom- 
etry for HD 133880 (Landstreet & Mathys 2000). This program has 
longitudinal and surface magnetic field strengths at four different 
phases as input, which are then used to iteratively search for the 
best-fitting i, p, B d , B q and fi ocl . For a more complete description, 
refer to Bailey et al. (2011). Measurements of ( B z ) from Land- 
street (1990), as well as the four new measurements from this work, 
provide good constraints on the observed variations. We also have 
important limits on (B) that appear to range from 10 to 20 kG near 
the positive and negative magnetic poles, respectively (see above). 
This information allowed us to derive a new field geometry: i 
55° ± 10°, 0 = 78° ± 10°, B d = -9600 ± 1000G, B q = -23200 ± 
1 000 G and B oct = 1900 ± 1000 G. Uncertainties were estimated by 
observing the deviations from the derived field by changing the input 
parameters of {B z ) and (£) within their uncertainties. This field cor- 



Figure 4. The top panel depicts the surface magnetic field modulus vari- 
ations as measured from Fe n (red upward facing triangles) and Cr ii (blue 
downward facing triangles), with the solid (black) curve representing the 
adopted best-fitting model. The bottom panel shows the (B z ) field variations 
observed for HD 133880 for our adopted magnetic field model. The red cir- 
cles are measurements from Landstreet (1990) using H /?, the blue triangles 
are LSD measurements from ESPaDOnS, the purple diamond HARPS and 
the solid (black) curve the adopted magnetic field model. The minimum of 
the (B z ) and {B) curves occur at phases 0.037 and 0.537, respectively. 

responds well with our calculated value of i = 63^ ± 1 8 and has the 
added advantage of agreeing with the alternate geometry proposed 
by Landstreet (1990) (see above). An illustration of our adopted ge- 
ometry is presented in Fig. 5. The bottom panel of Fig. 4 shows the 
model (B z ) curve overplotted with (B z ) measurements from Land- 
street (1990) and this paper. All the measurements agree well with 
the computed model and the measurements made using the metallic 
spectrum agree remarkably well with Landstreet (1990) who made 
measurements from Up. The field varies non-sinusoidally from 
about — 4.5kG at phase 0.037 to 2kG at phase 0.537. The field 
curve is flat and extended around the positive magnetic pole, with a 
rapid decline towards the negative magnetic pole. The quadrupole- 
to-dipole ratio, q, of the adopted model is 2.4, which clearly indi- 
cates the presence of a strong quadrupole component. We therefore 
concur with the conclusion of Landstreet (1990) that the field of 
HD 133880 departs in a significant way from a simple centred 
dipole. 

7 ABUNDANCE MODEL 

We utilize the fortran program zeeman (Landstreet 1988; Land- 
street et al. 1989; Wade et al. 2001) to perform abundance analysis 
of HD 1 33880. zeeman computes emergent spectra of all four Stokes 
parameters for stars with a permeating magnetic field. It assumes a 
low-order multipole geometry for the magnetic field based on the 
above derived values for i, P, £ d , B q and B ocX . Up to 10 phases of 
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Figure 5. The geometry adopted for HD 133880. The vertical axis is the rotation axis of the star. The angle between the line of sight and the rotation axis is 
i = 55". The angle between the rotation axis and the magnetic field axis is ft = 78 . The left-hand panel depicts an orientation where the positive magnetic pole 
is closest to the line of sight (</> = 0.537) and the right-hand panel is when the negative magnetic pole is nearer alignment with the line of sight = 0.037). 
The bands used for dividing the abundance distribution are shown as solid lines and are normal to the magnetic field axis. 


observed Stokes / spectra can be fitted simultaneously with the 
abundance of one element varied at a time. For each iteration, 
zeeman automatically adjusts vsin i and the radial velocity of the 
star to optimum values. Currently, up to six rings of uniform abun- 
dance on the surface of the star can be specified that have equal 
spans in magnetic colatitude, zeeman produces output for the abun- 
dance distribution of a given element at different magnetic latitudes. 
In effect, this model provides a rough ID map of the variation of 
abundance from one magnetic pole to the other. 

Initially, uniform mean abundances for each phase of HD 133880 
were found; however, we quickly discovered that a multi-ring model 
would provide a better fit. As explained by Bailey et al. (2011), 
zeeman was designed to quantify hemispherical abundance varia- 
tions for stars similar to 53 Cam (HD 65339; Landstreet 1988). 
HD 133880 exhibits large-scale abundance variations over the stel- 
lar surface and, as is evident from the longitudinal magnetic field 
curve, both magnetic hemispheres are observed. Our preliminary 
abundance analysis suggested that a maximum of three rings should 
be used to quantify the abundance variations because we saw abun- 
dance anomalies that varied between both magnetic poles and the 
magnetic equator. Nevertheless, we experimented with fewer and 
more rings on the surface of the star. We found that more than three 
rings did not significantly improve the quality of fits and fewer 
rings did not adequately model the observed variations in the Stokes 
/ spectra. 

For our analysis, T e( r was set to 13 000K and log g to 4.3, 
which are both consistent with the values determined in Section 4. 
The atomic data were taken from VALD (Piskunov et al. 1995; 
Ryabchikova et al. 1997; Kupka et al. 1999, 2000). The best-fitting 
value of o sin i was found to be 103 ± 10 km s -1 . Due to the strong 
blending in the spectral lines, identifying the best fit required vi- 
sual evaluation of the agreement between the computed and the 
observed Stokes / spectral fines in several windows of the order 
of 80-A wide. Since the strong magnetic field should suppress any 
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convective motions, which in any case are expected to be small at 
this r e fr, the microturbulence parameter was set to 0 km s“ l . 

7.1 Spectrum normalization 

For stars in which v sin i is greater than around 100 km s _1 (as is the 
case for HD 133880), the dominant source of error in abundance 
analysis is from continuum normalization. The reason for this is 
that in a spectral region that is of the order of 100- A wide, there 
may only exist a few continuum points from which to normalize the 
spectrum. Often it is difficult to determine which features should 
be placed near the continuum and which should not. To combat this 
problem, we performed continuum normalization of spectra of nor- 
mal, sharp-lined, A-type stars from both IheESPaDOnS andFEROS 
instruments to determine the highest degree polynomial required to 
normalize a given spectral window. All spectral windows were about 
100-A wide to maximize the number of continuum points present in 
the broad-line spectra. In all cases, no greater than two segments of 
a cubic spline were required to normalize the normal A-star spec- 
tra correctly. For a given instrument and spectral window, we also 
used two segments of a cubic spline to continuum fit the spectra for 
HD 133880. In this manner, we reduce the uncertainty introduced 
in the placement of the continuum and ensure that all spectra are 
normalized in the same way. Where necessary, we iteratively per- 
formed continuum fitting for spectral windows for which synthetic 
spectra from zeeman agreed poorly with our initial normalization. 
Hill (1995) discusses further challenges in performing abundance 
analysis on broad-lined stars. 


7.2 Choice of magnetic field model 

As described in Section 6, two magnetic field models adequately 
characterize the observed variations of the longitudinal magnetic 
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field, ( B z >: that of Landstreet (1990) and the one presented in this 
paper. Tc explore the effect of different magnetic field geometries 
on the derived abundances, we performed abundance analysis of 
HD 133880 for all elements using both geometries. The maxi- 
mum abundance differences between the two models were found 
to be 0.1 -0.2 dex. We therefore adopt the new magnetic field ge- 
ometry presented in this paper because it agrees well with the 
value of i derived from physical parameters and has the added 
constraints of the observed surface magnetic field modulus varia- 
tion, <B). We report all abundances using the new magnetic field 
geometry'. 

8 ABUNDANCE ANALYSIS 

A total of 12 spectra, well spread in phase, were available for HD 
133880, from the FEROS, ESPaDOnS and HARPS instruments. 
The HARPS spectrum was not used in the abundance analysis be- 
cause its phase corresponded closely to two other spectra with com- 
parable S/N and thus added no new relevant information. The large 
wavelength coverage of these instruments has allowed us to derive 
abundance distributions for O, Mg, Si, Ti, Cr, Fe and Pr, as well as 
upper limits for the abundances of He, Ni and Nd. The large v sin i 
presented challenges in deriving abundances due to the presence of 
significant blending in the spectra. Nevertheless, multiple relatively 
clean lines of the majority of elements studied were found with 
varying strengths and Lande factors, providing better constraints on 
the derived abundances. 

As discussed in the previous section, the present data set is suffi- 
cient to perform abundance analysis using a three-ring model with 
rings encompassing both magnetic poles and the magnetic equator 
with equal span in colatitude (60 ). Initially, six spectra were used 
(two ESPaDOnS and four FEROS) that were well spaced in phase 
to derive the three-ring abundance model. Once the mean abun- 
dance values were found for each ring, model fits were produced 
for all spectra. We have no observations for HD 133880 very close 
to the positive magnetic pole (phase 0.537). We do, however, have 
two spectra sufficiently near the positive pole, at phases 0.438 and 
0.640, adequate for providing a first approximation to the abundance 
distributions near the positive pole. 

The mean abundances derived for each element in each ring 
are tabulated in Table 6 along with solar abundances reported by 
Asplund et al. (2009). The quality of fits were tested (where pos- 
sible) using multiple spectral windows with derived uncertainties 
estimated from the observed window-to- window scatter of the abun- 
dance. For lines for which there was only one spectral line or region, 
we estimate the uncertainties to be about ±0.3 dex, which is the av- 
erage deviation from the derived abundance required to make the 
quality of fit unsatisfactory between the model and observed spec- 
tra upon visual examination. We also derived a set of abundances 


using T cf f = 12 000 K as reported by Netopil et al. (2008) for each 
of the elements listed in Table 6. We found that the abundances 
in each of the three rings changed by ±0.10 dex, or less, for the 
Fe-peak elements (Ti, Cr, Fe) and no more than ±0.1 5 dex for the 
rare earth elements (Pr and Nd), O, Mg and Si, which are within 
our quoted uncertainties. This suggests that the choice of 7^ does 
not significantly influence the derived abundances for HD 133880. 
To illustrate the quality of fits of the three-ring model, two 80- A 
windows for all 1 1 modelled spectra are shown in Figs 6 and 7. The 
quality and consistency of the model are discussed below for each 
element modelled. 


8.1 Helium 

Many possible lines of He i are considered when determining the 
abundance of helium: 447 1 , 5015, 5047, 5876 and 6678 A. However, 
none of these lines is unambiguously detected. Nevertheless, fitting 
these regions provides a useful upper limit to the abundance of 
helium. Our analysis confirms that HD 133880 is a He-wcak star 
with the abundance of helium at least a factor of 8 below the solar 
abundance. 


8.2 Oxy gen 

HD 133880 is classified as a Bp star with the Si A.4200 anomaly 
which suggests that the oxygen abundance should be near or less 
than the solar abundance (Sargent & Searle 1962; Roby & Lambert 
1990). Initially, we attempted to derive the oxygen abundance us- 
ing the 111 1-75 A triplet. These lines suffer from strong non-local 
thermodynamic equilibrium (non-LTE) effects and are highly sat- 
urated, making them unreliable for abundance determination. A 
more useful set of lines are those at 6155-56-58 A; the abundance 
distributions reported in Table 6 are from these lines. These lines 
are unambiguously detected, but are somewhat blended with Pr n 
near the redward wing. The O abundance varies from the magnetic 
poles to equator by a factor of about 2-3. At all phases, oxygen 
is found from these lines to be overabundant near the negative and 
positive magnetic poles by 0.5 and 0.3 dex, respectively, compared 
to the solar ratio. Near the magnetic equator, oxygen is comparable 
in abundance to the Sun, within uncertainties. 


8.3 Magnesium 

Only one clean line of sufficient strength to model was found for 
Mg ii at 448 1 A. This line shows some variability, implying an over- 
abundance of around 0.5 dex near both magnetic poles and a roughly 
solar abundance near the magnetic equator. 


Table 6. Abundance distribution of elements studied. Only upper limits of He, Ni and Nd are found. 

log(nx/«H) 

He O Mg Si Ti Cr Fe Ni Pr Nd 

Ring O' -60' (negative pole) <- 1.90 -2.78 -3.84 -2.79 -5.39 -4.55 -3.32 <-4.40 -6.52 -6.55 

Ring 60°-120 (magnetic equator) - —3.38 —4.37 —2.64 —5.62 —4.54 —3.44 - —6.72 —6.47 

Ring 120°- 180° (positive pole) -3.10 -3.91 -3.93 -6.75 -5.20 -4.11 -7.28 -8.22 

or - ±0.3 ±0.3 ±0.2 ±0.2 ±0.15 ±0.1 - ±0.2 ±0.3 

Solar abundance -1.07 -3.31 -4.40 -4.49 -7.05 -6.36 -4.50 -5.78 -11.28 -10.58 

# of lines modelled 2 1 1 4 3 3 3 1 3 1 
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Wavelength (A) 


Figure 6. Spectrum synthesis of the region 5000-5080 A using a three-ring model. The observed spectra are in black (solid) and the model fits are in red 
(dashed). From top to bottom, phases 0.054, 0.111 (both ESPaDOnS), 0. 1 39, 0.325, 0.438 (all FEROS), 0.640, 0.774 (both ESPaDOnS), 0.809 (FEROS), 0.878 
(ESPaDOnS), 0.917 and 0.990 (both FEROS). The line of eight is closest to the negative magnetic pole at phase 0.037 and the positive magnetic pole at phase 
0.537. Spectra are positioned at arbitrary positions along the vertical axis for display purposes. 


8.4 Silicon 

A total of four lines of Si n were found that were suitable for mod- 
elling: 5041 A, the doublet at 5055-56, 5955 and 5978 A. The first 
two lines of Sin were fitted simultaneously and the result tested 
by using the same derived abundance model in the latter two lines. 
Satisfactory fits were obtained at each phase for all lines of Si; how- 
ever, the doublet at 5055-56 A was systematically too strong in the 
model at all phases (see Fig. 6). Finding the best-fitting model for 
the longer wavelength lines did not rectify the problem of stronger 
spectral lines in the model and produced a similar abundance dis- 
tribution. We find large-scale abundance variations for Si on the 
surface of HD 133880 from pole to pole. At the negative magnetic 
pole and magnetic equator, Si is approximately 25 times more abun- 
dant than in the Sun. The lowest abundance is found at the positive 
magnetic pole where, however. Si is still about three times more 
abundant than in the Sun. 

Fig. 7 highlights three Sim lines at 4552, 4567 and 4574 A that 
are not modelled well using the abundance derived from the Si n 
lines discussed above (red-dashed line). An enhancement of the Si n 
abundance of about 1 dex is necessary to model well the Si ill lines 
(blue-dashed line). This discordance is very suggestive. The effects 
of non-LTE in the Si in lines are expected to be small at T cff less than 
about 15 000 K. Becker & Butler (1990) demonstrate that non-LTE 


effects in the Si m lines decrease towards lower temperatures. The 
most probable cause is strong stratification of Si in the atmosphere of 
HD 133880 with the abundance high in the atmosphere being much 
lower than near an optical depth (r c ) of about 1. However, the large 
v sin i (103 km s -1 ) makes normalization difficult (see Section 7.1) 
and studies on more hot Ap stars with low usin i values are necessary 
to determine how widespread this Si n/m discrepancy is, and how it 
may be explained. 


8.5 Titanium 

We have modelled the abundance of Ti using lines of Ti ii at 4563 
and 4571 A as well as at 4805 A. The final model was determined 
using the former two lines and tested using the latter line. TI exhibits 
drastic variations between magnetic poles of nearly 1.5 dex which 
can be seen in Fig. 7 where lines due to Ti ii are much weaker at 
phases 0.423 and 0.625 than near phase 0. Ti is most abundant at the 
negative magnetic pole where it is more than 40 times more abun- 
dant than in the Sun and transitions to slightly lower abundances 
at the magnetic equator (~25 times overabundant). At the positive 
magnetic pole, Ti is only about two times more abundant than the 
solar value. 
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Phase 

0.054 

0.111 

0.139 

0.325 

0.438 

0.640 

0.774 

0.809 

0.878 

0.917 

0.990 


4530 4540 4550 4560 4570 4580 4590 4600 

Wavelength (A) 

Figure 7. Same as Fig. 6 for the region 4530-4600 A. The green boxes highlight the three Si in lines at 4552, 4567 and 4574 A. The red-dashed lines indicate 
the model 5ts to abundances derived from the Si n lines that are presented in Table 6. The blue-dashed lines are model fits with about a l-dex enhancement to 
the Si abundances shown in Table 6. 

8.6 Chromium 

Several lines of Cr were found that were suitable for modelling, 
particularly the Cm lines at 4558, 4588 and 4592 A. The latter two 
of these three lines were relatively clean, while X4558 is blended 
with Fe n. All three lines were used to compute the three-ring model 
presented in Table 6. All Cr lines in Fig. 7 are fitted well at all phases 
with the model fit at 4592 A being slightly too strong at later phases. 

Large variations are observed from hemisphere to hemisphere. Cr 
is weakest at the positive pole, but its abundance is still 10 times 
greater than the solar abundance. The negative pole and magnetic 
equator have a mean abundance that is approximately 100 times 
greater than in the Sun. 

8.7 Iron 

A plethora of lines are available from which to derive the abundance 
distribution of Fe. The model was found using Fen lines at 4541, 

4555 and 4583 A. The lines were of sufficient strength for abundance 
analysis and only the line at 4555 A is blended with Cm. Similar 
to Cr, Fe is least abundant at the positive magnetic pole, where 
an abundance of about 2.5 times higher than the solar abundance 
is found, and most abundant near the negative magnetic pole and 
magnetic equator where a mean abundance approximately 10 times 
larger than in the Sun is observed. Most lines of Fe are modelled 
well, as seen in Fig. 7. Some weaker lines are systematically too 
weak in the model (see Fig. 6), which is most likely attributable to 
unrecognized blending. 

Z 2012 The Authors, MNRAS 423, 328-343 
Monthly Notices of the Royal Astronomical Society 2012 RAS 


8.8 Nickel 

There are a few possibly useful lines of Ni ii to model at 5058, 5059 
and 5064 A; however, none is unambiguously detected. Neverthe- 
less, the region (labelled in Fig. 6) proved useful for obtaining an 
upper limit to the abundance of Ni, which suggests that it is less 
than 10 times more abundant than in the Sun. 

8.9 Praseodymium 

Three lines of Pr were found that were useful for modelling. The 
final model was computed using Pr m and Prn at 6160 and 6161 A, 
respectively. These lines are blended with the redward wing of the 
O i triplet at 6155-56-58 A, but are of sufficient strength to obtain 
a useful abundance. The model fits were then tested using Ptiii at 
778 1 A. A mean abundance was found for the two rings covering the 
negative magnetic pole and magnetic equator that is about 1 0 5 times 
overabundant compared to the solar ratio. At the positive magnetic 
pole, Pr is least abundant but still 10 4 times more abundant than in 
the Sun. 

8.10 Neodymium 

A possible useful line of Nd m is at 6145 A; however, the S/N is such 
that it is not unambiguously detected at all phases, and rotational 
broadening causes this line to be strongly blended with Fe u at 6147 
and 6 1 49 A. We therefore turn to Nd hi at 47 1 1-1 2-1 4 A. These lines 
are relatively clean, providing useful abundance determinations that 


2 

1.8 

1.6 

1.4 

1.2 

1 

0.8 

0.6 

0.4 

0.2 

0 






l_ - __ V'l w " . I J I T T . . , 

u 











<n<n 


IX 


v/’V- 

£?sA 












^ ri 
<N <N . R CM U 
° 

__L Z. L 



JMZS* csr. 

<N m N m IN „ t. 

■- rT, C rT. .9 <N 


HD 133880 339 


reveal a strong variation in the mean abundance of Nd between 
magnetic hemispheres. For both the negative magnetic pole and 
magnetic equator Nd is > 10 4 times more abundant than in the Sun, 
whereas at the positive pole this rare earth element is about 200 
times the solar ratio. 

9 MAGNETOSPHERE 

Lim et al. (1996) demonstrated that the 3.5- and 6-cm radio flux 
and circular polarization of HD 133880 vary significantly and co- 
herently according to the ~0.877 d period. They reported that the 
emission shows broad peaks near the phases of the longitudinal 
field externa (corresponding, in Landstreet’s model, to the phases 
at which the poles of the dipole contribution to the field curve pass 
closest to the line of sight), and narrower peaks at the predicted 
phases of quadrupolc component pole passages (i.e. the quadrupole 
contribution to the field curve). They also reported that the strong 
extraordinary-mode circular polarization varied in sign with phase, 
in agreement with the sign of the longitudinal field. As described 
in Section 2, we have re-reduced the ATCA radio data, but there 
are no significant differences in the reduced data compared tc that 
described by Lim et al. (1996). Here we examine the radio emission 
in the context of a magnetospheric picture combining the dynamical 
concepts summarized by Owocki, Townsend & Ud-Doula (2006) 
and the radio magnetosphere model of Linsky, Drake & Bastian 
(1992). 

Using the stellar wind parameters determined using the CAK 
formalism (Castor, Abbott & Klein 1975: values of M = 
10“ 11 Mq yr -1 and = 750 km s -1 ), we compute the wind mag- 
netic confinement parameter tj* = B ^ R 2 /Mv 00 10 7 (neglecting 

clumping; ud-Doula & Owocki 2002), where B ^ is the equatorial 
surface strength of the magnetic dipole component. The rotation 
parameter W = v^lv orb — 0.2 (Ud-Doula, Owocki & Townsend 
2008). This places the Alfven radius at R^ = rj R ~ 60/?*. In 
contrast, the Kepler (or corotation) radius is located relatively close 
to the star, at about /?Ke P = W“ 2/3 /?* = 2.8/?*. This geometry, in 
which /?! cep <$C /?Aifi results in a large spatial volume in which mag- 
netospheric plasma is predicted to be maintained in rigid rotation 
with the star by the magnetic field while simultaneously being dy- 
namically supported against gravitational infall due to rapid rotation 
(Ud-Doula etal. 2008). 

In principle, this situation provides conditions suitable for signifi- 
cant accumulation of rigidly rotating stellar wind plasma (diagnosed 
from optical and ultraviolet (UV) line emission and variability; e.g. 
Petit et al. 2011), as well as a large region over which electrons can 
be accelerated (thus producing the non-thermal radio emission; e.g. 
Linsky et al. 1992). 

The phased radio light curves, illustrated in Fig. 8, exhibit a 
number of interesting characteristics. First, it is clear that the ra- 
dio measurements in Stokes I and Vfl vary significantly, and are 
coherently phased with the rotational period derived from the pho- 
tometric and magnetic measurements. Comparison of the 3.5- and 
6-cm polarization and light curves fully confirm the details of the 
variations. As described by Lim et al., the flux variation is com- 
plex, characterized by strong, broad maxima at phases 0.0 and 0.5 
(i.e. the extrema of the longitudinal field), and sharper, somewhat 
weaker secondary extrema at quadrature phases (i.e. 0.25 and 0.75). 
The polarization degree varies approximately sinusoidally, with ex- 
trema of hi 16 per cent and a mean of zero. The new rotational period 
derived in Section 3 brings the phases of the radio flux and polar- 
ization extrema into good agreement with the longitudinal field and 
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Figure 8. Lower frame: ATCA Stokes / flux variation at 6 cm (filled sym- 
bols) and 3 cm (open symbols). Upper frame: Stokes VU polarization at 
6 cm (filled symbols) and 3 cm (open symbols). 

photometric extrema, resolving the ~0.05 cycle offset pointed out 
by Lim et al. (1996). The approximate reflectional symmetry of the 
light curve in both I and VU (horizontally, relative to phase 0.0, as 
well as vertically for Stokes V7/, relative to the VU = 0) is in good 
qualitative agreement with the derived stellar geometry, in particular 
the large value of the magnetic obliquity p ~ 90°, implying that we 
view geometrically similar regions of both (magnetic) hemispheres 
during a stellar rotation (see also Fig. 5). 

Comparing the variations at 3.5 and 6 cm, we find that the ampli- 
tudes are nearly identical in both I and VII. The Stokes I variations 
may display some small differences in shape, in particular with the 
central (phase 0.5) peak at 3.5 cm being somewhat shallower and 
broader than at 6 cm. 

The presence of structures in the radio light curve attributable 
to the quadrupolar component of the field is very interesting. While 
the Alfven radius for the dipole component falls off as r) v J A y for 
the quadrupole the fall-off is as (Ud-Doula et al. 2008). The 
‘quadrupolar Alfvdn radius 1 for HD 133880, computed assuming 
the magnetic model derived in Section 6.2, is /? A if, q — 21/?*. This 
implies that the radio emission at both 3.5 and 6 cm forms signifi- 
cantly closer than 20/?* from the star’s surface, if the origin of the 
radio emission is charged particles in the outflowing stellar wind. 

The UV Civ A.X1548, 1550 and Siiv A.X1394, 1402 resonance 
doublets are both sensitive indicators of the stellar wind and its 
structuring and modulation by the magnetic field (e.g. Schnerr & 
Henrichs 2008). A search of the MAST archive indicates that HD 
133880 was unfortunately not observed by the IUE satellite, nor 
with any other instrument with coverage in this spectral region. On 
Lhe other hand, our optical spectra contain Ha as well as a number 
of Paschen series H lines. An examination of Ha reveals weak 
variability in the inner ±500 km s -1 of the profile. 

In Fig. 9 we illustrate the line profile variations that are observed 
in Ha. Fig. 10 shows the phased equivalent width of Ha, along 
with that measured from the Pa 8590 A line. A telluric line removal 
algorithm was applied to each profile before analysis to reduce the 
variations caused by such features. While particular care has been 
taken to keep the normalization of the ESPaDOnS and FEROS spec- 
tra as consistent as possible, it is difficult to judge whether the small 
fluctuations that are found in the far wings are real or an artefact of 
the normalization. At phases 0.0-0.35, we find that the line cores 
show excess absorption relative to the average profiles, reaching 
maximum absorption at phase 0.0. Between phases 0.35 and 0.8, 
the core region of the Ha profile shows emission relative to the av- 
erage profile. Unfortunately, our phase sampling is not sufficient to 
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Figure 9. Top: observed Her line profiles with the mean profile indicated 
as dashed red. Middle: shown are the continuum normalized spectra after 
subtracting the average of all the Ha line profiles, displayed in such a way 
that the profiles are plotted al a height that corresponds to the phase of the 
observation. Also included is a dotted line corresponding to a zero difference 
to highlight the line profile variations. Bottom: shown is the variance of the 
continuum normalized spectra from the mean profile. 

distinguish if these features propagate through the profile, as would 
be expected from rigidly rotating magnetically confined plasma. 

The Ha profile variability of HD 133880 might be a conse- 
quence of modifications of the atmospheric structure due to the 



Figure 10. Phased equivalent widths measurements of the core region of 
Ha (top) and the total Paschen line at 8590 A from the ESPaDOnS (black 
squares) and FEROS (red triangles) spectra. 

strong chemical non-uniformities we have determined to exist in 
its atmosphere. However, we know that the disc-integrated photo- 
spheric He abundance is constrained to be below 10 per cent that 
of the Sun, i.e. less than 1 per cent of hydrogen by number, at all 
phases. Therefore, even a significant variation of the He abundance 
could have no important effect on the atmospheric structure. 

For the He star a Cen (which has a huge variation of He abundance 
over its surface and strongly variable metallic lines; Bohlender, Rice 
& Hechler 2010), Ha varies only in the core of the line (over a region 
perhaps two or three times the width of a typical metallic photo- 
spheric line) and the variations quite closely mimic the variations 
observed in the star’s iron lines, i.e. the core of Ha gets stronger 
at the same time as Fe ii lines are strong. This may be what we are 
seeing in HD 133880 as well, since the surface field and metal lines 
are weakest when Ha is also weakest and most of the variations 
occur within ±200 km s -1 of line centre. Given the dominance of 
the quadrupolar field component, it is possible that the contribution 
of magnetospheric material to the overall variability of Ha may be 
reduced somewhat, since material will be trapped in a more complex 
circumstellar geometry: likely four or more ‘clouds’ rather than the 
two clouds seen in stars with dominant dipolar fields. For example, 
the Ha variability of HD 37776 (Grunhut et al., in preparation), a 
B-type star with a complex quadrupolar field, is quite modest com- 
pared to other stars with more dipolar fields (e.g. a On E; Oksala 
et al. 2011) despite the fact that it has similar physical properties to 
the other stars. Therefore, at present we simply conclude that the 
origin of the H line variations is unclear. 

10 DISCUSSION 

This paper presents our continuing work in the study of chemical 
abundance evolution with time of magnetic Ap/Bp stars in open 
clusters. The goal of this project is to model the magnetic field 
structure, to derive an estimate of the surface abundance distribu- 
tion of a number of elements and to provide a preliminary analysis 
to serve as a foundation for further more detailed modelling. In 
particular, more spectropolarimetric observations will allow for a 
more sophisticated interpretation of the surface abundance varia- 
tions and magnetic field structure via magnetic Doppler imaging 
(MDI; Piskunov & Kochukhov 2002). 

HD 133880 is a very young, rapidly rotating Bp star (v sin / ~ 
103 km s’ 1 ) that hosts a strong magnetic field. The magnetic field 
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variations were modelled using a colinear axisymmetric multipole 
expansion with i = 55° ± 10°, 0 = 78° ± 10°, = -9600 ± 

1000G, B q = -23200 ± 1000G and ft* = 1900 ± 1000G. We 
used Hipparcos photometric measurements together with our own 
{B z ) measurements, as well as photometric and radio emission data 
presented by Waelkens (1985) and Lim et al. (1996), respectively, 
to refine the period to P = 0.877 476 ± 0.000 009 d. The abundance 
distributions of He, O, Mg, Si, Ti, Cr, Fe, Ni, Fr and Nd were 
modelled using three coaxial rings around the two magnetic poles 
and the magnetic equator. The Fe-peak elements are overabundant 
compared to the solar ratios with strong variations on a global scale 
observed between hemispheres. 

We have 12 Stokes / spectra of HD 133880 spread well in phase 
that sample both magnetic poles and the magnetic equator. Our 
adopted magnetic field geometry, although an approximation to the 
actual field structure of HD 133880, and our three coaxial ringed 
abundance model provide the most detailed abundance analysis 
possible given the limitations of the model. Most elements studied 
appear to have a distinctly non-solar abundance in most rings. The 
abundance of Mg is slightly above solar at both magnetic poles 
but near solar abundance at the magnetic equator. Upper limits to 
the abundances of He and Ni were found. HD 133880 is clearly 
a He-weak star with an abundance at least a factor of 8 lower 
than the solar ratio. Ni is no more than 10 times overabundant 
compared to the Sun. All other elements studied (O, Si, Ti, Cr, 
Fe, Pr and Nd) are more abundant at the negative magnetic pole 
than the positive magnetic pole. Most interesting are the variations 
observed in O abundance. Most Ap/Bp stars show O abundances 
near or below the solar abundance ratios (Sargent & Searle 1962; 
Roby & Lambert 1990). However, a careful study of several O 
features reveals that oxygen is more abundant than solar at both 
magnetic poles (most notably at the negative magnetic pole where 
O is 0.7 dex overabundant), but a solar abundance near the magnetic 
equator. This behaviour is certainly uncommon in Bp stars, but is 
unambiguously determined in this analysis. 

In the context of the larger study of the time evolution of at- 
mospheric abundances of Ap/Bp stars, it is interesting to compare 
to work done by Bailey et al. (2011) on the Ap star HD 318107 
(NGC 6405 77). This star is also hot (7 cff ~ 1 1 800 K), young (log t 
~ 7.8 yrs), of comparable mass (2.95 ± 0.15 Mq) and host to a 
large complex magnetic field ({ B z ) varies from about 1 to 5kG). 
In contrast to HD 133880, it is a slow rotator (vsini ~ 7kms“ l ) 
and only one magnetic hemisphere is observed. The global prop- 
erties of atmospheric abundances is similar between both stars: the 
mean abundances for Fe-peak (Ti, Cr and Fe) and rare earth (Nd 
and Pr) elements are comparable within uncertainties and clearly 
overabundant compared to the solar ratios. There are, however, no- 
table differences in the abundances of O and Si between both stars. 
For HD 318107, the O abundance is more typical of Ap stars, be- 
ing about 0.3 dex lower than the solar value as compared to HD 
1 33880 where O appears to be overabundant. Si is about 1 dex more 
abundant in HD 133880 compared to HD 318107 which is under- 
standable in that the former has the Si A4200 peculiarity, whereas 
the latter does not. 

Folsom et al. (2007) studied the chemical abundances of stars 
in an older cluster NGC 6475 Qog t ~ 8.47 yrs) that have similar 
masses to HD 133880. Mean abundances for multiple elements 
in the Ap stars HD 162576 (NGC 6475 55; 3.0 ± 0.4 M 0 ), HD 
162725 (HR 6663; 3.3 ± 0.5 Mq) and HD 162305 (NGC 6475 14; 
2.7 ± 0.3 Mq) were derived. As compared to HD 133880 and HD 
318107, the abundances of Cr for all three stars are similar, being 
about 2 dex larger than in the Sun. However, the abundances of Ti 


and Fe appear to be, respectively, slightly lower and only 0.5 dex 
larger than the solar ratios. This is in contrast to HD 133880 and 
HD 318107 where Ti and Fe are, respectively, about 1 .5 and 1 dex 
overabundant compared to the Sun. A larger data set is required 
before any definitive conclusions can be made regarding the time 
evolution of abundances in Ap stars. 

IID 133880 was previously modelled by Landstreet (1990) who 
characterized the magnetic field variations of this star with mea- 
surements of the longitudinal magnetic field from H/L Unlike most 
magnetic Ap/Bp stars whose magnetic field structure is predomi- 
nantly dipolar, Landstreet (1990) found a magnetic field dominated 
by the quadrupole component. Despite the large v sin i of the star, we 
were able to further characterize the magnetic field of HD 133880. 
The higher S/N and spectral resolution of the ESPaDOnS spectra 
made it possible to compare spectral lines of the same element that 
have a large and small Lande z factor to characterize the surface 
magnetic field modulus variations, (5). By finding a model that 
characterizes the observed (B z ) and ( B ) variations, we are able to 
find the inclination, i, of the line of sight to the rotation axis as 
well as the tilt angle of the magnetic field axis, 0, to the rotation 
axis. This model is merely an approximation to the very complex 
field structure of HD 133880, but is sufficient to provide a mean- 
ingful first analysis to the atmospheric abundance distribution of 
elements. 

Fig. 4 shows clearly the qualitative nature of both the (B z ) and 
(£) variations. There is a rapid change in polarity in the {B z ) curve 
with a steep transition from the positive to negative hemisphere. 
The strong quadrupolar component of the magnetic field causes a 
distinct non-sinusoidal variation with the field structure flattening 
out near the positive magnetic pole. The fact that there is a clear 
pole reversal indicates that the sum of the inclination and magnetic 
field axis must be greater than 90 °J + 0> 90°, and that the line-of- 
sight goes well into both magnetic hemispheres. The parameters of 
the adopted magnetic field geometry are presented in Section 6.2. 
(ft) measurements using all metallic lines (see Table 2) agree well 
with the measurements by Landstreet (1990) using Hft suggesting 
that they sample the magnetic field in similar ways (see Fig. 4). 
Longitudinal field (B z ) measurements of spectral lines of only a 
single element provide clues to the inhomogeneous structure of 
the magnetic field and the surface distribution of elements. Using 
all four of our polarimetric observations, we have remeasured the 
longitudinal field ( B z ) using spectral lines of Si, Ti, Cr and Fe. The 
results are listed in Thble 7. 

It is clear from Table 7 that all four elements sample the field in 
different ways. The differences in measured field strengths can be 
attributed to the complex magnetic field (as reported by Landstreet 
1990, and discussed in this paper) and the inhomogeneous surface 
distribution of elements which is a common feature of Bp stars. 
Near the positive magnetic pole (phase 0.537) at phase 0.640, as 
well as phase 0.878, the measured (ft) values for all four elements 


Ihble 7. Measurements of {B z ) using LSD masks composed of lines of 
single elements (Si, Ti, Cr and Fe) from the polarized ESPaDOnS and 
HARPS spectra (cf. Table 2). 


Phase 

(B z ) (Si) 
(G) 

(B t ) (Ti) 
(G) 

(ft) (Cr) 

(G) 

(B z ) (Fe) 
(G) 

0.111 

—3419 ± 76 

-5724 ±535 

-4773 ± 229 

-3498 ± 50 

0.640 

1790 db 66 

2255 ± 304 

2653 ± 196 

1993 ±44 

0.787 

1020 ± 93 

-506 ± 308 

1026 db 208 

975 ± 78 

0.878 

-1776 ±63 

-2883 ±319 

-2409 ± 175 

-2010 ± 45 


2012 The Authors, MNRAS 423, 328-343 
Monthly Notices of the Royal Astronomical Society © 2012 RAS 



342 J. D. Bailey et al. 


differ significantly. In the negative magnetic hemisphere (near phase 
0.037) Si and Fe sample the field in similar ways, whereas Cr and 
Ti sample the magnetic field differently with (B z ) measurements 
approximately 1 and 2 kG stronger than Si and Fe. Along the sharp 
decline from the positive to negative magnetic pole (near phase 
0.787; see Fig. 4), it appears that only Ti differs greatly, showing 
a polarir reversal as compared to the other elements at the same 
phase. 

The measurements of the longitudinal field using different el- 
ements reveal that particular elements at different phases sample 
the field in very distinctive ways which suggests a truly complex 
magnetic field that cannot be accurately modelled by a simple ax- 
isymmetric magnetic field (such as the one adopted in this paper and 
by Landstreet 1990). It also strongly suggests an inhomogeneous 
distribution of various elements on the stellar surface that substan- 
tially affects the magnetic field measurements. We also compared 
the computed to the observed Stokes V profiles for the three po- 
larimetric ESPaDOnS spectra for multiple lines (Mg, Si, Ti, Cr and 
Fe). The agreement is poor, with the computed Stokes V profiles 
being much stronger than the observed, further suggesting that our 
adopted magnetic field model is only a rough approximation. This 
result is not unique. A recent paper by Kochukhov et al. (201 la) on 
HD 37776 (HIP 26742) reveals that the complex multipole model 
derived from the longitudinal field predicts Stokes V signatures 
that are ouch too strong compared to the observed Stokes V pro- 
files. In this case, (he multipole magnetic field model (similar in 
structure to the one used for HD 133880) predicts a large surface 
field, whereas the MDI results produce a much more modest surface 
field. This is also seen, although less severely, for a Ori E (Oksala 
et al. 2011). For this star, detailed analysis using MDI produced a 
magnetic field model that closely resembled the observed Stokes 
V profiles. The acquisition of new, high S/N polarimetric observa- 
tions of HD 133880 are necessary to characterize unambiguously 
the surface abundance variations and magnetic field structure using 
MDI. 

The very strong field and rapid rotation of HD 133880 produce 
interesting magnetospheric characteristics. In our preliminary anal- 
ysis, we were not able to conclude that the optical spectra con- 
tain detectable information about the magnetosphere. However, HD 
133880 is one of only a few Bp stars for which radio flux variations 
have been noted. There are small but distinct differences between 
the 3- and 6-cm radio emission, such as the systematically lower 
3-cm fluxes than 6-cm fluxes at phases 0.0 and 0.5. The 3-cm emis- 
sion suggests a ‘bump’ at about phase 0.6, and the peak at about 
phase 0.75 is narrower at 3 cm than at 6 cm. The fact that the radio 
spectrum is flat between 6 and 3 cm suggests that the emission is 
optically thick, and therefore we are observing the magnetosphere 
at two different radio surfaces rather than an integration throughout 
the entire volume. These subtle differences in the radio light curves 
at 6 and 3 cm then imply that these two surfaces cannot have ex- 
actly the same shapes and that, presumably, the magnetic field is 
not perfectly axi symmetric. HD 133880 is a very exciting target for 
more detailed magnetospheric modelling. 
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